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We solve the Sp(N) Heisenberg and SU(N) Hubbard-Heisenberg models on the anisotropic tri- 
angular lattice in the large-N limit. These two models may describe respectively the magnetic and 
' electronic properties of the family of layered organic materials ft-(BEDT-TTF)2X. The Heisenberg 

, model is also relevant to the frustrated antiferromagnet, CS2CUCI4. We find rich phase diagrams 

■ for each model. The Sp(N) antiferromagnet is shown to have five different phases as a function of 
the size of the spin and the degree of anisotropy of the triangular lattice. The effects of fluctuations 
at finite-N are also discussed. For parameters relevant to CS2CUCI4 the ground state either exhibits 

D ' incommensurate spin order, or is in a quantum disordered phase with deconfined spin-1/2 excita- 

»' 1 , tions and topological order. The SU(N) Hubbard-Heisenberg model exhibits an insulating dimer 

■ phase, an insulating box phase, a semi-metallic staggered flux phase (SFP), and a metallic uniform 
phase. The uniform and SFP phases exhibit a pseudogap. A metal-insulator transition occurs at 

, intermediate values of the interaction strength. 

S ' PACS numbers: 71.10.Hf, 71.10.Fd, 71.30.+h, 75.10.Jm 

o 

O \ I. INTRODUCTION 

The family of layered organic superconductors k-(BEDT-TTF)2X has attracted much experimental and theoretical 
*Z interest 

H|). There are many similarities to the high-T c cuprates, including unconventional metallic properties and 
competition between antiferromagnetism and superconductivity B. The materials have a rich phase diagram as a 
^1 ' function of pressure and temperature. At low pressures and temperatures there is an insulating antiferromagnetic 
£N) ■ ordered phase; as the temperature is increased a transition occurs to an insulating paramagnetic state. A first- 
order metal-insulator transition separates the paramagnetic insulating phase from a metallic phase; it is induced by 
increasing the pressure The metallic phase exhibits various temperature dependences which are different from 

that of conventional metals. For example, measurements of the magnetic susceptibility and NMR Knight shift are 
, consistent with a weak pseudogap in the density of states . 

The main purpose of this paper is to attempt to describe the magnetic ordering, the metal-insulator transition, 
and the unconventional metallic properties of these materials with two simplified models. We model the magnetic 
ordering in the insulating phase with the quantum Heisenberg antiferromagnet (HAF). We model the metallic phase as 
well as the metal-insulator transition with a hybrid Hubbard-Heisenberg model. To substitute for the lack of a small 
expansion parameter in cither model, we enlarge the symmetry group from the physical SU(2) = Sp(l) spin symmetry 
to Sp(N) (symplectic group) for the Heisenberg model (||]] an d to SU(N) for the Hubbard-Heisenberg model p0| , pl| . 
We then solve these models in the large-N limit and treat 1/N as our systematic expansion parameter. In section 
U we briefly summarize the physical SU(2) Heisenberg and Hubbard-Heisenberg model on the anisotropic triangular 



o 
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lattice. In section [II we review the large-N theory of the Sp(N) quantum Heisenberg model. Based on the large-N 
solution of this model, we present the magnetic phase diagram in the parameter space of quantum fluctuation n^/N 
(where rib is the number of bosons in the Schwinger boson representation of the spin) and the magnetic frustration 
J2I 'J\. We discuss the effects of finite-N fluctuations on the Sp(N) magnetic phases with short-range order (SRO). 
We also discuss how our results are relevant to understanding recent neutron scattering experiments on CS2CUCI4. 
In section |iy| we review the large-N theory of the SU(N) Hubbard-Heisenberg model. We present the phase diagram 
based on the large-N solution in the parameter space of the dimensionless ratio of the nearest-neighbor exchange to 
the hopping constant J\/t\ and the dimensionless anisotropy ratio J2HJ1 + J%)- Away from the two nested limits 
J\ = and J2 = we find a metal-insulator transition which occurs at finite critical value of J\jt\. We also find 
that the density of states in the metallic state is suppressed at low temperatures, in qualitative agreement with the 
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unconventional metallic properties seen in experiments. We conclude in section M with a brief review our results. 



II. HEISENBERG AND HUBBARD-HEISENBERG MODELS ON THE ANISOTROPIC TRIANGULAR 

LATTICE 



Based on a wide range of experimental results and quantum chemistry calculations of the Coulomb repulsion 
between two electrons on the BEDT-TTF molecules it was argued in reference || that the k-(BEDT-TTF)2X family 
are strongly correlated electron systems which can be described by a half-filled Hubbard model on the anisotropic 
triangular lattice. The Hubbard Hamiltonian is: 

h = -h J2 [4 a ^ + H - c -] - *2 Yl [ c i CTc J- + H - c i 

<ij> <<ij>> 

+ yE^-l) 2 (1) 

i 

Here c\ a is the electron destruction operator on site i and there is an implicit sum over pairs of raised and lowered 
spin indices a =\, \ . Matrix element t\ is the nearest- neighbor hopping amplitude, and ti is the next-nearest- neighbor 
hopping along only one of the two diagonals of the square lattice as shown in figure g. The sum < lj > runs over 
pairs of nearest-neighbor sites and << ij >> runs for next-nearest-neighbors. 




FIG. 1. The anisotropic triangular lattice with two types of bonds. Note that this can also be viewed as a square lattice 
with an additional next-nearest neighbor interaction along only one of the two diagonals. 



Physical insight can be attained by considering the Hubbard model for different values of the ratio U/t. In the 
limit of large U / 1 the Hubbard model at half-filling is insulating and the spin degrees of freedom are described by a 
spin- 1/2 Heisenberg antiferromagnet |f| : 

H = Ji £ £ • 5j + J 2 £ £.4 (2) 

<ij> <<ij>> 

where Si is the spin operator on site i, and the exchange couplings J\ = At\/U and J2 = 4i|/{7. Competition 
between J\ and J2 leads to magnetic frustration. Using parameters from quantum chemistry calculations fl^-|l5| it 
was estimated in reference || that J2/J1 ~ 0.3 to 1 for the k-(BEDT-TTF)2X family. Hence, magnetic frustration 
is important. The frustrated Heisenberg Hamiltonian interpolates between the square-lattice ( J 2 — 0) and the linear 
chain (Ji = 0). Much is known about these two limiting cases. Additional insight Jl6],[l7]] comes by considering 
different values of the ratio J2/J1. At J2 = 0, the square lattice limit, there is long-range Neel order with a magnetic 
moment of approximately 0.6^b; see reference Jl^ ]. If J2 is small but nonzero, the magnetic moment will be reduced 
by magnetic frustration. At J2/J1 around 0.5, quantum fluctuations combined with frustration should destroy the 
Neel ordered state. As J2/J1 is further increased, the system may exhibit spiral long range order p7y . At Jj = J2, 
the lattice is equivalent to the isotropic triangular lattice. Anderson suggested in 1973 that the ground state could 



be a spin liquid without long range order |18|. However, subsequent numerical work indicated that there is long- 
range-order with ordering vector q = (27r/3, 27r/3) |l1|. Finally, at J\ — the model reduces to decoupled Heisenberg 
spin- 1/2 chains which cannot sustain long range spin order [fl2j as per the quantum Mermin- Wagner theorem. For 
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Ji small but non-zero, the system consists of spin-1/2 chains weakly coupled in a zigzag fashion. The case of two 
such weakly coupled zigzag spin chains was studied by Okamoto and Nomura [^0| and by White and Affleck who 
showed that there is a spin gap in the spectrum A ~ e -constx J 2 /Ji^ an( j ^he ground state exhibits dimerization and 
incommensurate spiral correlations. Although our system consists of an infinite number of weakly coupled spin chains 
instead of two chains, we find similar behavior. 

As U/t decreases, charge fluctuations from electron hopping become significant. Competition between hopping 
and Coulomb repulsion leads to a transition from the insulator to a metal. We use the hybrid Hubbard-Heisenberg 
Hamiltonian pO| , pl| with independent parameters t, J, and U (where in general J ^ 4t 2 /U) to describe some aspects 
of the transition: 

H = E l-Wl"^ + H.c.) + MSi ■ § 3 - imnj)] 

<ij> 

+ E [-t2(4 a c ia + H.c.) + J 2 (S i -S i -^n i n i )] 

<<ij>> 

^E^- 1 ) 2 - (3) 

i 

Here n\ = c\° c\ a is the number of fermions at site i. The Hamiltonian reduces to the Hubbard model when J\ = J2 = 
and the Heisenberg model for tx,t2 — > 0. In the square lattice limit and at half-filling, perfect nesting drives the system 
to an antiferromagnetic insulator at arbitrarily small value of the interactions. As diagonal hopping ti is turned on, 
however, nesting of the Fermi surface is no longer perfect and the metal-insulator transition (MIT) occurs at non-zero 
critical interaction strength. 



III. SP(N) HEISENBERG MODEL 

We first focus on the spin degree of freedom, appropriate to the insulating phase of the layered organic materials, 
by solving the Heisenberg model in a large-N limit. Our choice of the large-N generalization of the physical SU(2) 
spin-1/2 problem is dictated by the desire to find an exactly solvable model which has both long-range ordered (LRO) 
and short-range ordered (SRO) phases. This leads us to symmetric (bosonic) SU(N) or Sp(N) generalizations. The 
former can only be applied to bipartite lattices. The latter approach has been applied to the antiferromagnetic 
Heisenberg model on the square lattice with first-, second-, and third-neighbor coupling (the J\ — J2 — J3 model) [0, 
the isotropic triangular lattice p2] j, and the kagome lattice ^2|. As the anisotropic triangular lattice is not bipartite, 
we must choose the Sp(N) generalization ||. 



A. Brief Review of the Approach 

To ascertain the likely phase diagram of the frustrated Heisenberg antifcrromagnet, we consider the Sp(N) symplectic 
group generalization of the physical SU(2)= Sp(l) antiferromagnet ||[22j- The model can be exactly solved in the 
N — > 00 limit. Both LRO and SRO can arise if we use the symmetric (bosonic) representations of Sp(N). We begin 
with the bosonic description of the SU(2) HAF, where it can be shown that apart from an additive constant, the 
Hamiltonian may be written in terms of spin singlet bond operators: 

H = -\ E Me a pb\ a bf)(e-< s b i7 b iS ) (4) 

ij 

where we have used the bosonic representation for spin operator 

Si = \b\ a 3*hf> , (5) 
where a =|, J, labels the two possible spins of each boson. The antisymmetric tensor, e a p is as usual defined to be: 
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We enforce the constraint rib = b\ b- la = 25 to fix the number of bosons, and hence the total spin, on each site. The 
SU(2) spin singlet bond creation operator e a pbl a b^ may now be generalized to the Sp(N)-invariant form J a pb\ a bY ■ 
Global Sp(N) rotations may be implemented with 2N x 27V unitary matrices U: 



b^Ub 
U^JU = J. 



(7) 



Here 



J 



V 



generalizes the e tensor to N > 1, and b- la with a 
of the Sp(N) HAF may then be written: 



(8) 



1 . . . 2N is the Sp(N) boson destruction operator. The Hamiltonian 



(9) 



where again the Greek indices run over 1 . . . 2N, and the constraint b\ a b- la = rib is imposed at every site of the lattice. 
We have rescaled Jy/2 — ► Jy/27V to make the Hamiltonian of O(N). For fixed J\j J% we have a two-parameter (rib, 
N) family of models with the ratio k = rib/N determining the strength of the quantum fluctuations. In the physical 
Sp(l) limit, k = 1 corresponds to spin-1/2. At large k (equivalent to the large-spin limit of the physical SU(2) model) 
quantum effects are small and the ground states break global Sp(N) spin-rotational symmetry. LRO then corresponds 
to Bosc condensation which we quantify by defining 



(10) 



birh,(T 



Here the spin-quantization axis has been fixed by introducing the paired-index notation a = (to, a) with m = 1 . . . N, 
a =|, I, and to — 2 ... N. The c- number spinors x a , when non-zero, quantify the Bose condensate fraction, and are 

. For sufficiently small k, however, quantum fluctuations overwhelm the tendency to 
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given by (b ima ) = 
order and there can be only magnetic SRO. 

Upon inserting equation [irj into equation |], decoupling the quartic boson terms within the corresponding functional 
integral by introducing complex-valued Hubbard-Stratonovich fields Qy directed along the lattice links, enforcing the 
constraint on the number of bosons on each site with Lagrange multipliers Ai, and finally integrating out the 6-fields, 
we obtain an effective action proportional to N. Thus, as N — > oo, the effective action may be replaced by its saddle- 
point value. At the saddle-point, the Qy, A,, and x- lr7 fields are expected to be time-independent, so the effective 
action can be put into correspondence with a suitable mean-field Hamiltonian. The Hamiltonian may be diagonalized 
by a Bogoliubov transformation, yielding a total ground state energy E MF given by p2[ 



E M f[QA] = £ / Jy_ |Q; . |2 _ ^ Q .. e ^ w + Hm 



N 



i k 



(11) 



where Wk are eigenenergies obtained from diagonalizing the mean-field Hamiltonian. Finding the ground state of the 
Sp(N) HAF now reduces to the problem of minimizing Emf with respect to the variables Qy and x- la , subject to the 
Lagrange- multiplier constraints 



dE MF [Q,\] 



. 



(12) 
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It is essential to note that the action possesses local U(l) gauge symmetry under local phase rotations by angle 
0i(r): 




,i9 I (T)+ie J (r) 




(13) 



This symmetry reminds us that the representation of spin operators in terms of the underlying bosons is redundant 
as the phase of each boson fields can be shifted by an arbitrary amount without affecting the spin degree of freedom. 
Two gauge invariant quantities of particular importance for our classification of the phases are |Qy|, and Aj. 

Extrema of Em f are found numerically with the simplex-annealing method |^3| . We work with a lattice of 40 x 40 
sites and check that this is sufficiently large to accurately represent the thermodynamic limit. Constraint equation |l2| 
is tricky, however, as Emf[M] is neither a minimum, nor a maximum with respect to the Ai directions at the saddle 
point. The problem is solved by decomposing Ai into its mean value A and the deviations from the mean, 5X{ = Ai — A. 
As A is gauge-invariant, we solve the constraint equation ^ separately for it by applying Newton's method. The 
resulting Emf[\ <5Ai] can then be maximized in the remaining <5Ai directions. Individual values of SXi are in general 
non-zero, but by definition it must be the case that ^> <5Ai = 0. 



To make progress in actually solving the model, we now make the assumption that spontaneous symmetry breaking, 
if it occurs, does not lead to a unit cell larger than 4 sites. Our choice of unit cell is shown in figure ||. 



FIG. 2. The 2 x 2-site unit cell for the Sp(N) Heisenberg antiferromagnet on the anisotropic triangular lattice. The com- 
plex-valued, directed, Qi to Q? fields live on the links of the square lattice, while Q\ to Q\ live on the diagonal links. Arrows 
denote the orientation of the Qy fields. 

The 2x2 unit cell requires 12 different complex-valued Q-fields (8 on the square links, and 4 along the diagonal 
links) and 4 different A fields and x spinors at each of the four sites. However, we have checked that at every saddle- 
point in the SRO region of the phase diagram (xi CT = 0) each of the 8 Q-fields on the horizontal and vertical links take 
the same value. Likewise, the 4 Q-fields on the diagonals are all equal, as well as all 4 A-fields. Thus the 2 x 2-site 
unit cell can be reduced to only a single site per unit cell as shown in figure |[ We expect this to also hold in the 
LRO phases, in accord with previous work on the Sp(N) model [p2~| . 



B. Phase Diagram of the Sp(N) Heisenberg antiferromagnet 
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FIG. 3. The one-site unit cell for the Sp(N) Heisenberg antiferromagnet on the anisotropic triangular lattice. 



The zero-temperature phase diagram is a function of two variables: J%j J\ and k. The various saddle points can be 
classified in several ways. Both SRO and LRO phases may be characterized in terms of an ordering wavevector q via the 
relation q = 2k m i n , where k m i n is the wavevector at which the bosonic spinon energy spectrum has a minimum. The 
spin structure factor 5(q) peaks at this wavevector 22 . LRO is signaled by non-zero spin condensate x- lai which we 
assume occurs at only one wavevector, fc m j„; that is, x^ = x exp(ifc m i„ ■ fyj and x^ — —ix exp(— ik m i n ■ fj). The phases 
may be further classified Q according to the particular value of q. There can be commensurate collinear ordering 
tendencies where the spins rotate with a period that is commensurate with the underlining lattice. Alternatively, 
there may be incommensurate coplanar ordering tendencies where the spins rotate in a two dimensional spin space 
with a period that is not commensurate with the lattice. 




FIG. 4. Zero-temperature phase diagram of the Sp(N) Heisenberg antiferromagnet at large-N. The strength of the quantum 
fluctuations is set by the parameter 1/k (vertical axis). Dimerization patterns induced by finite-N fluctuations are shown in 
the insets. 

The phase diagram is shown in figure 0. Note that the general shape of the phase diagram is qualitatively consistent 



with the finding of spin wave calculations [[17 16| that quantum fluctuations are largest for J2/J1 ~ 0.5 and large 
J2I J\- For large enough values of k, the ground state has magnetic LRO. As the magnetic LRO phases break Sp(N) 
symmetry, there are gapless Goldstone spin-wave modes. As a check on the calculation, we note that in the k — > 00 
limit there is a transition between the Neel ordered and incommensurate (q, q) LRO phase at J2/J1 = 0.5 in agreement 
with the classical large-spin limit. At smaller values of k there are quantum disordered phases that preserve global 
Sp(N) symmetry. In the N — > 00 limit these are rather featureless spin liquids with gapped, free, spin-I/2 bosonic 
excitations (spinons). Finite-N fluctuations, however, induce qualitative changes to the commensurate SRO phases 
(see below). In the limiting case of the nearest- neighbor square lattice, J2 = 0, we reproduce the previously obtained 
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result |g| that Neel order arises for k > 0.4. In the opposite limit of decoupled one dimensional spin chains, J\ = 0, 
the ground state is in a disordered phase at all values of k. There are five phases in all, three commensurate and two 
incommensurate, as detailed in the following two subsections. 



1. Commensurate phases 

At small to moderate Jij J\ there are two phases with Q\ x = Qi y 7^ 0, and Q2 = 0. The eigenspectrum Wk has its 
minimum at k = ±(7r/2, 7r/2), with the implication that the spin-spin correlation function peaks at q = (tt,tt). Neel 
LRO with x- la 7^ appears when k is sufficiently large. The boundary between LRO and SRO phases is independent 
of J2/J1 except at one end of the boundary, but this is expected to be an artifact of the large-N limit ]23). Finite-N 
corrections should bend this horizontal phase boundary. 

At large J2/J1 and small k the ground state is a disordered state characterized by Q\ x = Q\ y = and Q2 7^ 0. The 
chains decouple from one another exactly in the large-N limit, but at any finite-N the chains will be coupled by the 
fluctuations about the saddle point. Also, x\„ = as it must, by the Mermin- Wagner theorem. The Sp(N) solution 
does not properly describe the physics of completely decoupled spin-1/2 chains. All spin excitations are gapped, and 
there is dimerization at finite-N (see below). This behavior is in marked contrast to the gapless, undimerized ground 
state of the physical SU(2) spin-1/2 nearest-neighbor Heisenberg chain jjjj. 



2. Incommensurate phases 

At intermediate values of J2/J1, there are two incommensurate phases with Q\ x = Q\ y 7^ 0, and Q2 7^ 0. The 
ordering wavevector q = (q, q) with q varying continuously from 7r to 7r/2, a sign of helical spin order in a given plane. 
The inverse k~ x of the critical k separating LRO from SRO peaks near the isotropic triangular point J2/J1 = 1, 
where it agrees in value to that reported in a prior study of the triangular lattice p2] ] , and then decreases with a long 
tail as J2/J1 increases. As the one-dimensional limit of decoupled chains is approached, n c — > 00. Again this accords 
with the Mermin- Wagner theorem. 

All the phase transitions are continuous except for the transitions between the (jr, tt) LRO phase and the (q, q) SRO 
phase which is first order. Fluctuations at finite-N, however, modify the mean-field results |^). The modifications are 
only quantitative for the LRO phases, and for the incommensurate SRO phase. In particular, instantons in the U(l) 
gauge field have little effect in the incommensurate SRO phase which is characterized by non-zero Q%. The Q2 fields 
carry charge ±2, so when they acquire a non-zero expectation value, this is equivalent to the condensation of a charge- 
2 Higgs field. Fradkin and Shenker showed some time ago that a Higgs condensate in 2+1 dimensions quenches 
the confining U(l) gauge force. Singly-charged spinons are therefore deconfined, instead the instantons which carry 
magnetic flux are confined, and no dimerization is induced || . The relevant non-linear sigma model which describes 
the transition from an ordered incommensurate phase to a quantum disordered phase has been studied by Chubukov, 
Sachdev, and Senthil @. 

In the case of the commensurate (71", tt) SRO and the decoupled chain phases, instantons do alter the states quali- 
tatively. The Berry's phases associated with the instantons lead to columnar spin-Peierls order [^4],[| (equivalent to 
dimer order) as indicated in figure H. The dimerization pattern induced in the decoupled chain phase is similar to 
that found by White and Affleck |21jj for a pair of chains with zigzag coupling. Furthermore, spinon excitations are 
confined into pairs by the U(l) gauge force. 

Note that the deconfined spinons in the (q, q) SRO phase are qualitatively different than the spinons found in the 
limit of completely decoupled chains, J\ = 0. The (q,q) SRO phase exhibits true 2-dimensional fractionalization, 
in contrast to the decoupled one-dimensional chains. The spinons are massive, again in contrast to those found in 
one-dimensional half-odd-integer spin chains. The transition from the dimerized chain phase at small J1/J2, which 
has confined spinons, to the deconfined incommensurate phase at larger J\/ J% is described by a 2+1 dimensional Z2 
gauge theory |27]]. In fact the (q,q) SRO phase is similar to a resonating valence bond (RVB) state recently found 
on the isotropic triangular lattice quantum dimer model ps| . The phase has "topological" order; consequently when 
the lattice is placed on a torus (that is, when periodic boundary conditions are imposed), the ground state becomes 
four-fold degenerate in the thermodynamic limit p9,BQ|. 
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C. Physical Spin- 1/2 Limit 



It is interesting to examine in more detail the physical spin- 1/2 limit corresponding to k — 1. In figure]^ we plot the 
ordering wavevector q as a function of the ratio J 2 /( J\ + J 2 ). Note that quantum fluctuations cause (i) the Neel phase 
to be stable for larger values of J2/J1 than classically, and (ii) the wave vector associated with the incommensurate 
phases deviates from the classical value. Similar behavior was also found in studies based on a series expansion fl3l| , 
and slave bosons including fluctuations about the saddle point u2j. 
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FIG. 5. Ordering wavevector q (in units of 7r/a) of the large-N Sp(N) Heisenberg antiferromagnet at k = 1, which corresponds 
to spin-1/2 in the physical Sp(l) limit. The solid line is the classical ordering wavevector p7| . 

Commensurate q = ir Neel order persists up to J 2 /(Ji + J2) = 0.369, and this is followed by the incommensurate 
LRO phase for 0.369 < J 2 /(Ji + <h) < 0.880. Then a tiny sliver of the incommensurate SRO phase arises for 
0.880 < J 2 /(Ji + J 2 ) < 0.886. Finally there is the decoupled chain phase for 0.886 < J 2 /(Ji + J 2 ) < 1. A strikingly 
similar phase diagram has been obtained by the series expansion method 0. A comparison between the Sp(N) and 
series expansion results is shown in figure |6[ Both methods suggest that there exists a narrow SRO region between 
Neel and incommensurate ordered phases, though in the Sp(N) case the region does not extend all the way up to 
k = 1. Possibly finite- N corrections to this large-N result could change the phase diagram quantitatively such that 
the (q,q) SRO phase persists up to n = 1 in agreement with the series expansion results. Another difference is that 
the Sp(N) result shows no dimerization in this narrow SRO region (due to the above-mentioned Higgs mechanism), 
while the series expansion indicates possible dimer order. 

(7C, 7C) (q, q) (q, q) Dimer 

LRO LRO SRO SRO J, 

SP(N) I 1 H 1 

0.0 0.369 J 1.0 J x + J 2 

0.880 0.886 
no dimerization 



(^ I) Dimer (<J> q) 

LRO SRO LRO J 2 

Series Expansion | 1 — | 1 1 

0.0 0.411 0.473 0.8 1.0 J x + J 2 
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FIG. 6. Comparison of Sp(N) phase diagram at k — 1 with the results of a series-expansion |B]J 



Similar results have also been obtained in a weak-coupling renormalization group (RG) treatment of the Hubbard 
model on the anisotropic triangular lattice p3| . For the special case of a pure square lattice (with next-nearest- 
neighbor hopping <2 = 0) at half-filling, the antiferromagnetic (AF) couplings diverge much faster during the RG 
transformations than couplings in the BCS sector, indicating a tendency towards magnetic LRO. As £2 is turned on, 
BCS and AF instabilities begin to compete. For sufficiently large <2, a crossover occurs to a d x 2_ y 2 BCS instability, 
suggesting that the system is now in a magnetic SRO state. Further increasing ti to reach the isotropic triangular 
lattice (t\ — t%) there are indications that long-range AF order re-enters. Finally, for £2 3> t\ LRO tendencies are 
again eliminated, this time by the strong one-dimensional fluctuations, ft is remarkable that the same sequence of 
ordering and disordering tendencies - LRO to SRO to LRO to SRO - occurs in all three approximate solutions. 



D. Application To Materials 

The region of the phase diagram intermediate between the square lattice and the isotropic triangular lattice is 
relevant to the insulating phase of the organic ^-(BEDT-TTF^X materials. The expected range in the spin exchange 
interaction j| is J2/J1 ~ 0.3 to 1. Depending on the precise ratio J2/J1, our phase diagram indicates that these 
materials could be in the Neel ordered phase, the (q, q) LRO phase, or possibly the paramagnetic (q, q) SRO phase 
(see above). In fact antiferromagnetic ordering with a magnetic moment of 0.4 to 1 x fis per dimer has been seen 
in the splitting of proton NMR lines in the k-(BEDT-TTF) 2 Cu[N(CN) 2 ]C1 compound @. It is conceivable that a 
quantum phase transition from the Neel ordered phase to the paramagnetic (<?, q) SRO phase or to the (q, q) LRO 
phase can be induced by changing the anion X. 

Coldea et al. |54| have recently performed a comprehensive neutron scattering study of CS2CUCI4. They suggest 
that this material is described by the spin-1/2 version of our model with J2/J1 = 2.5 and J2 = 0.37 meV. The 
measured incommensurability with respect to Neel order, 7r — q, is reduced below the classical value by a factor 0.47, 
consistent with series results J3lj (our Sp(N) solution shows a smaller reduction). Maps of the excitation spectra 
show that the observed dispersion is renormalized upward in energy by a factor of 1.67, which can be compared to 
theoretical values of 1.18 for the square lattice and 1.57 for decoupled chains. Furthermore, the dynamical structure 
factor S(q,ui) does not exhibit well defined peaks corresponding to well defined spin-1 magnon excitations. Instead 
there is a continuum of excitations similar to those expected and seen in completely decoupled spin-1/2 Heisenberg 
chains. In the case of a chain these excitations can be interpreted as deconfined gapless spin-1/2 spinons. 

In the relevant parameter regime, J2/J1 — 2.5, the large-N Sp(N) phase diagram predicts that the ground state 
is spin-ordered with an incommensurate wavevector (q,q). However, as noted above, finite-N corrections could move 
the phase boundaries so that the physical spin-1/2 limit is in fact described by the (q,q) SRO phase. In this phase 
there is a non-zero gap to the lowest lying excitations [which occur at wave vector (k m i n ,k m i n ) = (g/2,q/2)] rather 
than gapless excitations. At k = 0.56 and J2/J1 = 2.5, in the large-N limit, the system is in the (q,q) SRO phase 
of figure |] and the spinon gap is approximately 0.05 Ji. This is much smaller than the resolution of the experiment 
[see Figure 2(a) in reference j34j], in which excitation energies have only been measured down to about 0.5 Ji. The 
corresponding spinon dispersion is shown in figure 0. Again we stress that the spinons in the incommensurate SRO 
phase are qualitatively different than those which arise in a one-dimensional chain. 
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FIG. 7. Spinon dispersion along the (1, 1) direction in the deconfined (q, q) SRO phase for k — 0.56, Jl = 0.148 meV, and 
J2 = 0.37 meV; thus, J2/J1 = 2.5. The lowest energy excitations occur at an incommensurate wavevector of k m 0.26 x n 
where there is a non-zero energy gap of approximately 0.02 meV. 

In the deconfined (q, q) SRO phase, there are no true spinwave excitations, as spin rotational symmetry is unbroken. 
Nevertheless, when the gap to create a spinon is small, the spinwave description remains useful. For example, in 
neutron scattering experiments, spinons are created in pairs, as each spinon carries spin-1/2. So a spinwave may be 
viewed as an excitation composed of two spinons, though of course this spinwave does not exhibit the sharp spectral 
features of a true Goldstone mode. At large-N the spinons do not interact; finite-N fluctuations will lead to corrections 
in the combined energy and to a finite lifetime. The minimum energy of such a spinwave of momentum q is given, in 
the large-N limit, by: 

E q - = Mm{uj(q/2 + p/2) + oj(q/2 - p/2)} (14) 

where the minimization is with respect to all possible relative momenta p. The resulting spinwave dispersion is plotted 
in figure || alongside the classical result, scaled to the same value of the spin. The Sp(N) calculation shows a rather 
large upward renormalization in the energy scale compared to the classical calculation; this is a result of the quantum 
fluctuations which are retained in the large-N limit. 
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FIG. 8. Minimum energy of two spinons in the (q, q) SRO phase (solid line) compared to that of classical spinwaves (dashed 
line) for the same parameters as in figure (?]. The dispersion, in the case of the Sp(N) theory, is obtained from equation 
The classical spinwave dispersion is obtained from reference jl7| where it was scaled to spin-1/2. These energies have been 
multiplied here by a factor of k 2 « 0.31 to account for the reduced magnitude of the spin. The large upward renormalization 
in the Sp(N) excitation energies compared to the classical energies is due to quantum fluctuations. 

In the incommensurate (q, q) LRO phase (k = 1 and N — > oo) the spinon spectrum has gapless excitations, as 
shown in figure ^|. Apart from the absence of the small gap, the spinwave dispersion in the ordered phase is similar to 
the minimum energy spectrum in the disordered phase. Again there is an upward renormalization of the energy scale 
as shown in figure [L0|, though the ratio is relatively smaller than in the more quantum k = 0.56 case. The size of the 
renormalization is in good agreement with that seen in the CS2CUCI4 experiment It is important to note that 
finite-N gauge fluctuations bind spinons in the LRO phase into true spinwave excitations, with corresponding sharp 
spectral features. In contrast, as noted above, spectral weight is smeared out in the SRO phase. A large spread of 
spectral weight is seen in the neutron scattering experiments |mJ . But as the incommensurate SRO and LRO states 
are separated by a continuous phase transition (see figure |j), in the vicinity of the phase boundary it is difficult to 
distinguish the two types of excitation spectra. The spin moment in the LRO phase is small there, as is the gap in 
the SRO phase. Further experiments on CS2CUCI4 may be needed to determine which of the two phases is actually 
realized in the material. 



0.8 




(k,k) x 71 

FIG. 9. Spinon dispersion along the (1,1) direction in the incommensurate (q,q) LRO phase. As in figure 0, Jl = 0.148 
meV and J2 = 0.37 meV, but now k = 1. Gapless excitations occur at an incommensurate wavevector of k ~ 0.27 x tt. 
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FIG. 10. The same as figure g except in the incommensurate (q, q) LRO phase with the same parameters as in figure ^: 
Jl = 0.148 meV and J2 = 0.37 meV, and k = 1. 

IV. SU(N) HUBBARD-HEISENBERG MODEL 

We now turn our attention to the charge degrees of freedom by studying the hybrid Hubbard-Heisenberg model. 
This model should provide a reasonable description of the layered organic materials because the Hubbard interaction 
is comparable in size to the hopping matrix elements, U w t. Thus the stringent no-double-occupancy constraint of 
the popular t — J model should be relaxed. In the large-N limit it is also better to work with the hybrid Hubbard- 
Hciscnberg model than with the pure Hubbard model because the crucial spin-exchange processes are retained in 
the large-N limit [^lj. In the Hubbard model these are only of order 1/N and therefore vanish in the mean- field 
description. 

As there are now both charge and spin degrees of freedom, it is no longer possible to employ purely bosonic variables, 
in contrast to the previous section. Instead we use antisymmetric representations of the group SU(N) as the large-N 
generalization of the physical SU(2) system. As shown below, this generalization precludes the possibility of describing 
magnetic LRO or superconducting phases, at least in the large-N limit. However, as the same representation is placed 
on each lattice site, this large-N generalization works equally well for bipartite and non-bipartite lattices. It has been 
applied to the spin- 1/2 Heisenberg antiferromagnet on the kagome lattice J36|| . 

A. Brief Review of the Approach 

The Hubbard-Heisenberg Hamiltonian on the anisotropic triangular lattice is specified by equation ^. The gener- 
alized SU(N) version is obtained ]Q] by simply letting the spin index a in equation |^ run from 1 to N (where N is 
even) : 

H = [-*i( c i CTc j CT + H.c.) - ^{c\ a c ia cfc ip + -n s rij)] 

<ij> 

+ Yl [ _ *2(c| CT Cjv + H.c.) - '^{c^c^cfcip + ^nirij)} 

<<ij>> 

+ ^E( c N<,-^/2) 2 (15) 

i 

where all spin indices are summed over. Here we have also rescaled the interaction strengths Jj/2 — ► Jj/iV and 
[7/2 — > U /N to make each of the terms in the Hamiltonian of order N. At half-filling, the only case we consider here, 
a further simplification occurs as the term J\$n\n- } is simply a constant. There is no possibility of phase separation 
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into hole-rich and hole-poor regions, nor can stripes form p7|, as the system is at half-filling. We drop this constant 
term in the following analysis. 

We could also include a biquadratic spin-spin interaction of the form J(Si ■ Sj) 2 with J > 0. In the physical 
SU(2) limit this term does nothing except renormalize the strength of the usual bilinear spin-spin interaction. But 
for N > 2 it suppresses dimerization |ll[] , as the concentration of singlet correlations on isolated bonds is particularly 
costly when the biquadratic term is included. Thus there exists a family of large-N theories parameterized by the 
dimensionlcss ratio J/J, each of which has the same physical SU(2) limit. In this paper, however, we set J = as 
we find that the phase most likely to describe the metallic regime of the organic superconductors is a uniform phase 
with no dimerization which is stable even at J = 0. 

After passing to the functional-integral formulation in terms of Grassman fields, the quartic interactions are de- 
coupled by a Hubbard-Stratonovich transformation which introduces real-valued auxiliary fields <p on each site and 
complex- valued \ fields directed along each bond: 

0. = i £ (crCia _7V/2), Xij = ^ c - CjCT . (16) 

Clearly (</>i) is proportional to the local charge density relative to half- filling (corresponding to N/2 fermions at each 
site) and \ij may be viewed as an effective hopping amplitude for the fermions between site i and j. The effective 
action in terms of these auxiliary fields, which may be viewed as order parameters, can now be obtained by integrating 
out the fermions. We note that, unlike the bosonic formulation of the pure antiferromagnet, here there is no possibility 
of magnetic LRO as the order parameters x and 4> are global SU(N) invariants, and of course there is no possibility of 
Bose condensates in the fermionic antisymmetric representation of SU(N). Superconductivity likewise is not possible 
in the large-N limit because the order parameters are invariant under global U(l) charge symmetry rotations. In 
the Hciscnbcrg limit t — > the action is also invariant under local U(l) gauge transformations as long as the \~ and 
^-fields transform as gauge fields: Xij( r ) ~* e^ ei ( T ^ Si ( T ^Xiji T )i <M T ) — * <M T ) — d9i(r)/dT. For the more general case 
of Hubbard-Heisenberg model, this local U(l) gauge symmetry breaks to only global U(l) gauge symmetry reflecting 
the conservation of total charge. 

Since S e ff has an overall factor of N, the saddle-point approximation is exact at N — > oo. We expect cj> and x to 
be time-independent at the saddle point, so S e ff can be written in terms of the free energy of fermions moving in a 
static order-parameter background: 

S eff[<t>,x] = (3F[(f>,x;fA , 

F\&)Cj4 _ \ - IXijl 2 IXijl 2 V-r 1 , 2 _ * , , 

<ij> <<ij>> i 

" ^E ln { 1 + ex P[-^-^)]} ■ (17) 
' k 

Here the Wk are the eigenenergies of the mean- field Hamiltonian Hmf- 

<ij> <<ij>> i 

In the zero-temperature f3 — > oo limit the fermionic contribution to the free energy reduces to a sum over the occupied 
energy eigenvalues. The saddle point solution is found by minimizing the free energy with respect to x fields, and 
maximizing it with respect to the <p fields. We carry out the minimization numerically via the simplex- annealing 
method pj| on lattices with up to 40 x 40 sites. 

After ascertaining the zero-temperature phase diagram we then study the effects of non-zero temperature. As the 
temperature is raised, (3 — > and the last term in equation [T^ approaches ks In 2 per site reflecting the fact that each 
site is half-occupied. The entropy then dominates the free energy, terms linear in Xij i n equation ^ disappear, and 
the free energy is minimized by setting Xij = 0- Thus as the temperature is raised, antiferromagnetic spin correlations 
are weakened and then eliminated altogether. 



B. Zero- Temperature Phase Diagram of the SU(N) Hubbard-Heisenberg Model 

We again choose a 2 x 2 unit cell, in anticipation that translational symmetry can be broken at the saddle points. 
The 2x2 unit cell requires 12 different complex x-fields and 4 different real 0-fields as shown in figure [llj At half-filling 



13 



all fa = 0. As expected, there is no site-centered charge density wave, and the phase diagram does not depend on the 
size of the Hubbard interaction U (so long as it is repulsive) because fluctuations in the on-site occupancy, which are 
0(y/N) N/2 are suppressed in the large-N limit. Therefore the saddle-point solutions may be classified solely in 
terms of the remaining order parameter, the x-fields. For the special case t\ = t-2 = it is important to classify the 
phases in a gauge-invariant way because there are many gauge-equivalent saddle-points. In this limit there are two 
important gauge- invariant quantities: 

(i) The magnitude |xij| 2 which is proportional to the spin-spin correlation function (Si • S$). Modulations in \x\ signal 
the presence of a bond-centered dimerization. 

(ii) The plaquette operator II = Xi2X23X34X4ij where 1, 2, 3, and 4 are sites on the corners of a unit plaquette. By 
identifying the phase of x as a spatial gauge field it is clear that the plaquette operator is gauge- invariant, and its 
phase measures the amount of magnetic flux penetrating the plaquette 1 1 1 . Different saddle points are therefore 
gauge equivalent if the plaquette operator has the same expectation value, even though the x-fields may be different. 
In the Heisenberg limit, the flux always equals to or tt (mod 2ir), so a gauge can always be found such that all the 
X-fields are purely real. 




FIG. 11. The 2 x 2-site unit cell of the SU(N) Hubbard-Heisenberg model on the anisotropic triangular lattice. Arrows 
denote the orientation of the complex- valued Xij fields. 



Away from the pure Heisenberg limit t\ = t-2 = wc further classify the saddle-point solutions in terms of whether or 
not they break time-reversal symmetry (T). Finally, as there are four independent parameters (ti, £2, Ji, and J2) the 
phase diagram lives in a three-dimensional space of their dimensionless ratios. To reduce this to a more manageable 
two-dimensional section, we assume that J1/J2 = (ii/^) 2 ; then by varying J\ and the two hopping matrix elements 
we explore a two-dimensional space. The resulting phase diagram is shown in figure |l2|. We summarize the phases 
which appear in the diagram in the following subsections. 
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FIG. 12. The zero temperature phase diagram of the SU(N) Hubbard-Heisenberg model on the anisotropic triangular lattice. 
The metallic uniform region is indicated by dark shading while light shading demarks the semi-metallic staggered flux phase. 



1. One-Dimensional Dimer Phase 



This phase exists in the region J2 > Ji > t\ and it exhibits spin-Pcicrls (= dimer) order. All X2-fields are negative 
real numbers with \\ = xh x\ = x\i an d IX2I > Ixll- Also, all xi-fields are either small negative real numbers or 
zero with \\ = Xi = Xi = Xi < 0) x\ — Xi — Xi — Xi = an d IXil ^ Ixll- See figure [l3| for a sketch. The system 
breaks up into nearly decoupled dimerized spin chains. This phase breaks preserves T-symmetry as all the \~ fields 
are real. It is insulating as there is a large gap in the energy spectrum at the Fermi energy. The phase is very similar 
to the decoupled chain phase of the insulating Sp(N) model; in fact the dimerization pattern is identical to one of 
two such possible patterns in the Sp(N) model (see figure f|) and echos the pattern found by White and Affleck for 
the two-chain zigzag model ^J. In the extreme one-dimensional limit of J\ = t\ = the SU(N) solution, like the 
Sp(N) solution, fails jll| to accurately described the physics of decoupled chains, for at half-filling, the physical one- 
dimensional SU(2) Hubbard model is neither dimerized nor spin-gapped. The inclusion of the biquadratic interaction 
suppresses dimerization |38| and yields a state qualitatively similar to the exact solution of the physical system, but 
for simplicity we do not consider such a term here. 



FIG. 13. One-dimensional dimer phase. Dark lines indicate bonds with strong spin-spin correlations. These are the dimers. 



2. Box (also called "Plaquette") Phase 



This phase is also insulating and consists of isolated plaquettes with enhanced spin-spin correlations [|39|,[30| . See 
figure || for a sketch. All xi-fields are complex with \ X {\ = \xl\ > Ixfl = lx!l > IXil = Ixfl > Ixfl = IxiT The X2 
fields are small, | X2 1 = Ixll = Ixll = Ixll ^ \xl\- The phase 9 of the plaquette product X1X1X1X1 is neither nor 
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7r in general. The box phase breaks T-symmetry when t\ =^ 0. As time-reversal symmetry is broken, there are real 
orbital currents circulating around the plaquettes |hi| as shown in the figure. Apart from T-breaking, this phase is 
rather similar to the (n, n) SRO phase of the Sp(N) model as it is a commensurate SRO phase with a large spin gap. 




FIG. 14. Box (or "plaquette") phase with circulating currents. 



3. Staggered Flux Phase (SFP) 

All %1-fields are equal, with an imaginary component, and the X2-ficlds are equal, real, and much smaller in 
magnitude than the xi fields. The phase of the plaquette operator differs in general from or 7r; see figure ^ for a 
sketch. The staggered flux phase breaks T-symmetry. Like the box phase, there are real orbital currents circulating 
around the plaquettes pOfl , in an alternating antiferromagnetic pattern. The SFP is semi-metallic as the density 
of states (DOS) is small, and in fact vanishes linearly at the Fermi energy at J 2 = 0. Apart from T-breaking this 
phase is rather similar to the (7r,7r) LRO phase of the Sp(N) antiferromagnet, as the spins show quasi-long-range 
order with power-law decay in the spin-spin correlation function. In the limiting case of a pure Heisenberg AF on 
the nearest-neighbor square lattice, ti = t 2 = J 2 = 0, gauge fluctuations at sufficiently small-N are expected to drive 
the SFP (which can be stabilized against the box phase by the addition of the biquadratic interaction) into a (n, w) 
Neel-ordered state [|l)- 
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FIG. 15. Staggered flux phase (SFP) with circulating currents. 



4- Uniform Phase 

All x-fields are negative real numbers and Xi = Xi = ' ' ' = Xit X2 = X2 — X2 — X2 with xi x\ m general. 
See figure |l6| for a sketch. This phase preserves T-symmetry and since all x _ fields are real, they simply renormalize 
hopping parameters t\, The uniform phase is therefore a metallic Fermi liquid. Spin-spin correlations in the 
uniform phase decay as an inverse power law of the separation, with an incommensurate wavevector, so the uniform 
phase behaves similarly to that in the (q,q) LRO phase of the Sp(N) model (see below). 
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FIG. 16. Uniform phase with no broken symmetries. The uniform phase is a Fermi liquid. 

C. Global Similarities Between SU(N) and Sp(N) Phase Diagrams 

There are some global similarities between the SU(N) phase diagram of the Hubbard-Heisenberg model (figure |lj) 
and the Sp(N) phase diagram of the insulating Heisenberg antiferromagnet (figure ||). We have already pointed out 
similarities between the four phases of the SU(N) model and the decoupled chain, (tt, n) SRO, (it, it) LRO, and (q, q) 
LRO phases of the Sp(N) model. Apparently the dimensionless parameter J\jt\ on the vertical axis of the SU(N) 
phase diagram (figure |l|) is the analog of the quantum parameter 1/k of the Sp(N) phase diagram. The reason 
why these two dimensionless parameter play similar roles can be understood by considering the limit of the pure 
insulating antiferromagnet, corresponding to £1 — ► and £2 — > 0. In these limits, the SU(N) model can only be in 
the purely insulating box phase, or in the dimerized phase. The spins are always quantum disordered, and behave 
like the extreme quantum k — ► limit of the Sp(N) model. In the opposite limit £i — > oo and t% — ► oo the spin-spin 
correlation function decays more slowly, as an inverse power law instead of exponentially. This is as close to LRO as 
is possible in the large-N limit of the SU(N) model. Roughly then it corresponds to the ordered classical n — > oo limit 
of the Sp(N) model. 

D. Observable Properties of the SU(N) Hubbard-Heisenberg Model 

We now comment on the possible relevance of our large-N solution to the electronic and magnetic properties of the 
organic superconductors. In reference || it is estimated that J2/J1 ~ 0.3 to 1 and J\jt\ ~ 0.5 to 2. Figure [l2] then 
implies that the ground state is the uniform phase, which as noted above, is a rather ordinary Fermi liquid with no 
broken symmetries. 

We note that in the one-dimensional limit (t\ = J\ — 0) and in the square lattice limit (£2 = J2 — 0), nesting of the 
Fermi surface is perfect, and the system is driven into an insulating phase no matter how large the hopping amplitudes. 
Away from these two extreme limits, however, there is a non-trivial metal-insulator transition line. This accords with 
expectations because as the Fermi surface is not perfectly nested, the metallic state is only eliminated at a nonzero 
value of J\/t\. Experimentally it is found that upon increasing pressure, a SIT transition from the antiferromagnetic 
insulator to a superconductor is seen in the k-(BEDT-TTF)2X family of materials (3^^^]. This transition can be 
understood in terms of the SU(N) phase diagram as follows. As pressure increases, the effective hopping amplitudes 
£1 and £2 also increase because Coulomb blocking is less effective ||. The bandwidth broadens, but the spin exchange 
couplings Ji and J2 remain nearly unchanged as these are determined by the bare, not the effective, hopping matrix 
elements. Thus the ratio Ji/£i is reduced at high pressure, and the many-body correlations weaken. From the phase 
diagram (figure |l^) it is apparent that the system can be driven from an insulating state into a metallic state, which 
presumably superconducts at sufficiently low temperature. Thus the transition to a conducting state at high pressure 
can be seen as being due to bandwidth broadening, as in the Brinkman-Rice picture of the MIT. 
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FIG. 17. Density of states per unit volume at the Fermi energy as a function of temperature for three points, indicated by 
the inset and specified in table M, inside the uniform phase. 



TABLE I. Parameters (in meV) for the three points in figure hjl 



Symbol In Figure 17 


h 


h 


Ji 


J 2 


Ji/h 


J 2 /(Ji + J 2 ) 


star 


6 


5.4 


24 


19.2 


4 
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We plot the DOS as a function of temperature in figure |17| for three points inside the uniform phase identified in 
table |. The DOS decreases as the temperature decreases from room temperature down to absolute zero. Clearly 
the pseudogap is more prominent near the boundary between the conducting uniform phase and the semimetallic 
SFP and insulating box phases. The drop in the DOS could qualitatively explain the 30% depression seen in the 
uniform susceptibility (3f| (NMR experiments find a reduction of about 50% as the temperature decreases from 
100 °K to 10 °K HQ). The behavior may be understood as follows. As the temperature decreases, antiferromagnetic 
spin-spin correlations develop. These correlations are signaled in the mean-field theory by the link variables Xij which 
become non-zero at low enough temperature. Electrons on neighboring sites then tend to have opposite spin, reducing 
Pauli-blocking, and making their kinetic energy more negative. The band widens and the density of states drops. 

V. DISCUSSION 

We have solved the bosonic Sp(N) Heisenberg and fermionic SU(N) Hubbard-Heisenberg models on the anisotropic 
triangular lattice in the large-N limit. The bosonic Sp(N) representation of the Heisenberg model is useful for 
describing magnetic ordering transitions. It therefore may be an appropriate description of the insulating phase of the 
layered organic superconductors k-(BEDT-TTF)2X and of the insulating CS2CUCI4 compound. The fermionic SU(N) 
Hubbard-Heisenberg model provides a complementary description of the charge sector, in particular the physics of 
metal-insulator transitions and the unconventional metallic phases. Systematic expansions in powers of l/N about 
the large-N limit are possible for either model. 

For the Sp(N) model we found five phases: (i) three commensurate phases, the (71", tt) LRO and SRO phases, and 
the decoupled chain phase and (ii) two incommensurate phases, the (q,q) LRO and SRO phases. Passing from the 
square lattice (J2 = 0) to the one-dimensional limit (Ji = 0) at k = 1 which corresponds to spin-1/2 in the physical 
Sp(l) limit first there is the Neel ordered phase, the incommensurate (q,q) LRO phase, the incommensurate (q,q) 
SRO phase, and finally a phase consisting of decoupled chains. These phases are similar to those obtained from 
a series expansion method |]3~i|| and a weak coupling renormalization group technique |33| . The effects of finite-N 
fluctuations on the saddle-point solutions were also discussed. The observed dispersion of spin excitations in the 
CS2CUCI4 material p4[ ] can be described either as spinwaves in the incommensurate (q,q) LRO phase, or in terms of 
pairs of spinons in the deconfined (q, q) SRO phase. 

The zero-temperature phase diagram of the SU(N) Hubbard-Heisenberg model in the large-N limit has a 1-D dimer 
phase, a box (or plaquette) phase, a staggered-flux phase, and a uniform phase. In the extreme 1-D and square 
lattice limits, the ground state of the half-filled Hubbard-Heisenberg model is always insulating because the nesting 
of the Fermi surface is perfect. But away from these two extreme limits there is a metal-insulator transition. For 
parameters appropriate to the k-(BEDT-TTF)2X class of materials we find that the conducting regime is described 
by our uniform phase, which is a rather conventional Fermi liquid with no broken symmetries. In this phase we find 
that the density of states at Fermi level decreases at low temperatures, due to the development of antiferromagnetic 
correlations. This could explain the depression seen in the uniform susceptibility of the organic superconducting 
materials at low temperatures. It agrees qualitatively with experiments which suggest the existence of a pseudo-gap. 

We have used two models instead of one because the two large-N theories have complementary advantages and 
disadvantages. The bosonic Sp(N) approach is not suitable for describing electronic properties, as there are no 
fermions. The fermionic SU(N) approach, on the other hand, is not useful for describing magnetic ordering, as the 
order parameters are all SU(N) singlets. Neither mean-field approximation can describe the superconducting phase. 
The fermionic Sp(N) Hubbard-Heisenberg model does support superconducting states, but no non-superconducting 
metallic phases |p7| . Whether or not a single model can be constructed which is exact in a large-N limit and yet 
encompasses all the relevant phases remains an open problem. 
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